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B^.n ^liloli taiea 'lAto aoequat-- th-e.-aet'SCs^jitfij'o 

^f!!."** S° ^^'^ deterainatloa of the torsional. 

Indloftt. tHat the aerodynanic epan- ef fect -nay iaor,aa> •' 
^^^^ ry*"*"! reloom^B fou^d ^ neai^a of tha s.ction- 
roru« theory as much as 17 to 40 percent. 

It le found that the magaltud© of the effect in- 
lllW "ilu le&ree of ctiffa«ss taper and d«- 

creaees vlth increasing degree of chord taper;f^ ;,:v.'-f \ ? 

By a 83?ght sxtenslon of the pree?nt jn©^hc4 in 
pOBBiMe to ;:.ia.ly-j4 the el.astic dei-ora..-, ',loi,a of wings 
and the resul^xant .lift distributions, tofo.-e- torsltfftai ■ 
divergence occura . 



Thl» paper deals vith 'the linitlng ease of the l>©nd- 
ing-toraion flutter problea vhich occur* when tfe«.^f2:w^trer' 
frequency ha« the rafue zero. Thi« aspect of. the pro^iea 
Iia«J>een formulated and dealt vith, ^i^-ai-pi'^ltii^^i^^'-- - 
i*?%i%5oc***'S*^ ia«fratility. by H. Eel<,en«r (r*,f^r8»r.« 
tL: *r*&ti^«at'.t. wt»li i5e't.i5.«*4Ri*r - 

• wotk on the general flutter prohlea by TheoAartjen (?cf~ 

! expr«»8l6as. fixv'tfeo Y^isvaat air f-V ■> . -.^ 

\ aeotlon of the vthr vhlcii eo77«ff9«n& to 
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The purpose of the present paper is to inveetigate 
the effect of this eimplifying assumption by giving a 
procedure for the analysis of the t or sional-divergence 
problem which takes into account the aerodynaaio span . 
effect. The development 9 are based on the theory of tor- 
sion of straight rods and on lifting-line theory for the 
spanwise distribution of lift. 

A rapidly convergent process of iteration is devised ^ 
for the solution of the equations of the two theories for 
an elastically twisted wing. The method is applied to 
some typical examples and it is found that for a wing 
with an aspect ratio of about six the aerodynamic span 
effect modifies the torsi onal -divergence velocity obtained 
with the assumption of air forces of the two-dimensional 
theory by 17 to 40 percent, depending on the elastic and 
plan-form characteristics of the wing. 

This invest igation', conducted at the Massachusetts 
Institute of Technology, was sponsored by, and conducted 
with financial assistance from, the Sat ional Advisory 
Comai.ttQe for Aeronautics, 

- v-i ; .: .o " ' ' 

SYMBOLS 



b wing span 
c wing chord 
Cj^ roo^ chord (c* » c/o^) 

angle of attack before elastic deformation 
* angle of twist due to elastic deformation 
I section lift (per unit span) f^f^' ^^q&^^t^ 

p,., > density: Of. air.. - - 
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P auxiliai^y lift ;fuDptlon .^ccj/mc^, ^ ^<^^)^J^^r ^-^ 

e diftance be tween ,o_enter.^of -prea sure .and^ 

ei^ value of e at root (e* » e/e^) 

G .. modulus of ri.gidi,ty ilL.. . ' . . .. 

Cfl torsional rigidity of wing flection^**/ 



Ir value of I at root (I* = I/Ir)' ' ' / -^^-v 

S ■< '-Tprp^eot'ed wing "area" ^ ^-^ . ' \i ^.^t ' . 

v-l^a ;\ dimendlonlesfl conetant : -(mCR/4b ) . . r .: " v? / ' 

■ " J /mp eiD Vc^ A* " 
P torelonal-dive^T-gence parameteri / =^ — ^ V| ' : 

y -»'r^^tfJ|; r^jjauwiee cdordlnat 6 measured from wing root in 

C^i*/i^^^aMdMLIlary:^ 'defined in B:}\xatiQnn ' (SSYi^^i^^f:^'^-^^ 

-vtv^?,^;:jV;Hg7), and (^rX . ^. ^ v.- > v- .: .- ^^^^^ 

H^fi^^^ tf^s'^a'nf ^ def i^Vift taper 'ch'aracter'i «t4^# 
At'tC arbitrary,, constants 

^q^^ yWj^ilB^ coef fi ci <int •^ofJ^Simpaon^ fi ^tfle 
?f corresponding to section-force theory (ae subecript) 




•If a wing with ^an . initial angle-of-attack distribu-- 
tion ao(y) is BUbJ«<^ted 1 6 a"aif t dietribut: on l(y) 

^p#r ustit span >nd:?tli^^r Insult ant ^Ir .forest Idy^ aBsooi- 
ated with a npaiwiee^v^lematit' of the wl^agvacte at a dla- 
tauoe from the elaetic axis of the wing, a chan^c« of 
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angle of attack J takes place due to the torsional flex 
i\)llity of the wing. (See fig. 6.) The relationship be- 
tween the additional angle of attack and the lift per 
unit span is given approximately by the differential 
equation . , 

e(y) Ky) (l) 



where y is a spanwif^e coordinate measured from the wing 
root in units of half span h/2, e is the distance 
between the center of pressure and the elastic axia» and 
GI is the torsional rigidity. It is known that this ap- 
proximate equation neglects the effect of the spanwise 
variation of twist on the gtreesea and deformations of 
the wing, and that in a core accurate theory equation (l) 
would be replaced by a fourth-order differential equation 
for t). Since, however, all . calculations of the diver- 
gence velocity by means of the section-force theory of 
which the authors have knowledge are based on equation 
(1) and the main purpose of this paper is the estimation 
of the aerodynamic span effect, if is thought that satis- 
factory results may be obtained if the present calcula- 
tions . also are based on this equation. It may be stated 
that there are no essential difficulties in. extending the 
work of this paper in the direction of refined procedures 
for the determination of the elastic def ormati one • 

An additional relationship is afforded by the lifti;ig 

line integral equation ^ 

.(,) » .CW{^^' (a, * » I - -^f- ^} 

where c is the chord of the wing, p is the density of 
air, V is the velocity of flight ,. and m, la a profile 
constant. The notation f' i» used In equation' (2) to 
indicate that the Cauchy principal value of the integral 
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If an auxiliary lift function r(y) is defined ae • 
r(y) . . ^•c;(y) (3) 

where ot(y) is the conventional a ect 1 on-llf t co'ef f f - ^ ' ^ 
cient, equations (l) and (2) can he written in ths fOiMB-^ 

A. \i*{y)^]^ e*(y) F(y) « 0 (4) 
.dy L dy J 



c 



:kl4.t / ^ i:^. , aoCy) * «(y) (5)' 
*(y ) TT dtl y-H 



vrhe'r#^^'^--|i'^-'ir a. diaeneionresa parameter 



and ^7 relationship 
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In these equations c^ represents the root chord while 
o*(y) is the ratio of the chord to the root chordi 

c * ( y ) = ' V ^ 

• . .. ' T i^^i Cr> . . r^v: .:i^i^ 

and analogous definitione apply to ej, e*, Ij. and 

* . .. -r : .'':;v^-'; r . . ^- ' • ■ ^ ' i ^ ? ''-^'.'r^ t''2'"r i?'<?^'?-lL?'|:^.^^-- 
" ^ i^'j fv ,tW; F^-^^^ twisting\ atttba>>r^^»^ii 
(y 0):, , a<i4rn^t.aJi*^ tvjtatinc «ioak«nte are -applied, the ril- 
^ ouxk4j»^ry:r e on dl t lo n e f o^;* .„t t e f Txn jp 1 1 ^ a?# i- " 
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i{0) = 0 



(8) 



= 0 



while vaniehdng. of the lift at 
"boundary oondittone 



the wing tips leads to the 



F(±l) = 0 



(9) 



The determination of continuous* functions 7 and 
•3 satisfying the s inultaneoufl equations (4) and (5), 
together with the "boundary conditions of equations (8) 
and (9), constitutes the basic prohlem of torsional di- 
vergence. It should he rpaarked that unless ^ ia an 
odd function of y the three boundary conditions of 
equation (8) do not»in general., permi t a regular solution 
for The physical explanation of this occurrence lies 

in the fact that unless the wing loading is antisymmet- 
rical with respect to the wing root the restraint offered 
hy the' fuselage , in this formulation of the problem, is 
equivalent to a concentrated twisting moment, so that 

the function I — must have ^ discontinuity at the 

dy 

root (y = 0), "Also, since it is necessary that the 
function F(y) have a continuoue derivative at inte- 
rior points of the span in order that the left-hand side 
of Equation (5) he continuous, it follows that the func- 
tion 




which represents the so-called "induced angle of attack*... 
must have a di scont iuuous first derivative at the wing y;.'* 
root . \: : ' • 

^ • ■ ' [ . . . • ✓ . . 

Equations (4j and (5) can be combined into & single 
integro-dlf f erent lal equation by^ th6 ellminatloh'-'bfj^he 
function d. Thus, if equation (4> 1 s integrated tWitfe 
and the boundatr^' oanditicn« of equatlbn' (6)' 1 
an expression for i cau be determined in. the fora 
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. A •v.^.v'i" ^'^ - - aV'^* f * - * s-.> '^ht;:^ vi l^.V's.- t-^iS5;* , 

vh'6re the' fUnetionT" 0- depends upon the funCtl^on I',. 

the IntVoduotloiv of :equatioV^ iix.to ea\Maition\.(5i;^iV^ik 

• • • ' , ■ ■ - :\ -;(^ 

If .the function is not identically eero, equa- 

tion (ll), together with the "boundary conditions of equa-\ 
tloa (9);^ poBeeaseSvin general, a-untque^aolUtion ' ''T(^y:y^^^V[ 
It la known, however, that there exists an infinite set 
of critical values of the paraueter 3 for which no 
sqXutlqin to-equation: (ll) exists. In , addition ^.-a^ ■^'-.y'^-'^^ 
approaches one of these critical valuer , the magnitude i?^ " 
the corresponding functions ^(y) and * (y ) . increaeea 
without limit. Since 3 is proportional to the velocity 
of flight, the critical values of P correspond to crit- 
ical velocities at which a very large (theoretically in- 
finite) twisting force is experienced hy the ^wing.^ Thue 
an accurate determination of the smallest critical value 
of 9 Is deslrahle for purposes of structural wing de- 
el^gzi;. The value of V .corresponding to the smallest 
oritlcal valuesof 3 le designated aa the torsional- 
dlvar^ane« velocity^ / V 

According to the theory of integro-differentlal equa- 
tions the critical values of for which no solution to 
equation (ll) exists are Ide&Vlcal with the values of 9 
for which the homageneous equation, with Identically. 

rWa, poseesees a solution. That is , ^ a critical value of 
& oorreeponds to such a critlcaI:;Lveloclty of flight that 
an initially untwisted airfoil KflL ^[O) may become de- 

formed. In the present llnearlted theory the magnitude 
' ofv tl^v.dff'ljiCtlonvln: thlft^ ;Oae«t le^ ofc ttntdetarmlii^d Jfta^^ 
tudA *eoi«t!r*. iOft ti^^^^ nhm 

val^isrnof the 'parame^ g ■ l^^jt^i^:^^^ :.{n . . 
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" 8eotlon-f ore© theory," which diarttgards the aerodynamic 
effect of finite span "by neglecting the integral reDre- 
senting the induced angle of attack: in eauatlon (5) (ref- 
erence !)• In what follow5?,a hrief treatment of the 
section-force procedure i? firot given, after -which a 
method of successive approximation is presented for the. . 
determination of the t or si onal-diyergence 'flight veloc- . 
Ity and of the form of the corresponding deflection and . 
lift curves, according to the lifting-line theory '.of* 
equations (4) and (5), Since the magnitude of the 
t orslonal-divergence velocity is independent of the ini- 
tial angle-of-attack distrihution, it will be assumed' 
that the wing is initially at a tero angle of attack 
(a^(y) = 0). 

SOLUTION OF TH J PROBLEM jLCCOEDIHO TO SSCTION-FORCS THSORY 



If the integral in equation (6) is neglected and an 
Initial zero angle of attaclc is assumed, equations (4) 
and (5) "become \ 

ri*(y) ^1 ^ P%-(y)F(y) 0-'^^f-;'^?.^\Cl2) 
dy I dyj ' ' ' ■ "-^ ' 



These equations, together with the boundary conditions 

* (0 ) = 0 C 
. ; f «; ^ ^ a^^'i y^tK^ .^^^ 

are taken' a« • t he baei s: of? the at aly s 1 • - ofti >th«r ixir'oVle*" of - - 
torsional divergehce accordin'^ to: the' secti^n-f o^rW the^-. 
ory. The boundary conditions of equation (9) are not 
prescribed in this theory. 



to 



Vhich together with tieTKaio^feil^^ "boundary oonditions 
of equation (14) l8 etifficlent' Vo determine the critical 
values of the parameter P and the corresponding critical 
deflection modes,' (See reference !•) 

^'^^ A Class of S^pllclt Solution* 

The integration of equation (15) in closed form is 
possible, in particular, in,.oaa^,a. vheav^J^.^ .pho,r* 
section stiffness vary according to the laws 

c'(y) » e*(y) » (l - ay) • I 

I^(y) « (1 - ay) 

^ nUiere .'^j: /Vat^c-a^^ ^l^r>^tt positive re"5'ijstAn*«. ^ 
a It- a ^po^'sitiYe constant "less than unity. With .the ^,.iuV 
ititutlon' .; , > V \l' '^fi-. 

>^t|rOn (15) heooii^ee S v»ry V 

-U«v^:^;: . . . . '^z'-" • V;.'- :=■ 

the solution ie knoirn to ha oVv thtii :for^ 



19) 



where X , . is the ^e;i.eral^ ^Pea sai fuaction.of order V 

and ^' •• " ' ' ^ ' " M 

and . , « c'r *v '.s ' ) . ? in v: 



2u 

. 6 =. 



(31) 



- tn the special cases" Vhen , tf. :ro t :■■ . 

Ya = 2 (Yi + 1) (22) 

the solution can 'be expressed ia the form 

♦^y)'«W"'^'- ^ ^"yA^sln (X log yi") B coa- CX Ibg yiiV :-i435) 
Where A and B are arbitrary constants and 




The solution ia. evaa.\^a^ted' eicj)li<;itly In the follow- 
"Ing four cai'e«-: cr ^ ' • v ( : - ; •• --^ "* v [ . 

1. Uniform chor d» uniform etiffneef (V^ - '^a * Q) #** 
In this case equation (15) becomee "iij j^a^ 

- s; , (25) 

and the general continuous solution having eontinu one ^ 
derirativee exc^fpt at y s^-O-^d^f^f the form 



11- 



where the flrat term haa a discontinuous derivative at the 
root . The ■boundary aondltions of • .equatlott; (I4)- Xea^VJfvJ>? 
that ' '■• . 

and . 

cob' P « 0 

Squation (28) indicates 'that equatione (25) and (14) 
poeseea solutions only when 

_ i2n 1) tt- ';:rA^ii^^-; V^' (39^, 

Where n is an integer,^ The enallest of these values, 
3 3 tt/2, then corresponds to the t orsional-divergenoe 
velooity, and the corresponding deflection ao^^© .^^ J^^6. 

• ^#«i> -^(7> • i: Bin Plyl B sin Py " 

where A and B are arbitrary constants. 

Xouation ,(30) shove that this mode »dy have both 
eymaetrioal and ant i symmetrical component a, sd' that ac- 
cording to this theory the two halves of the wing- deflect 
independently of each other. Moreover, according to the 
Beotlt>ix-tforca .th-eory ^he symmetrical and antleyM^tr^X: : 
deflection modes correspond to the same critical fligftVv- 
velocity. Tor this reason it will he ooiivenient in this 
section to c oneider Conlys Cone^half of ^he, wing. The solu- 
tion of the problem for the first deflection mode then 

6 » 



A-V ^ ... § • • 

a. Uplfor« chogd. (tuadrat lolly decrc atinff <ti^fa««i 

•ld»;«d. tht een«r&l tolutioa 'of^ •o.ufctioii Ci6) i. d^^tlTnara 
frea equfttton (83) Xn the foT» (t«« «ito r«f«r*»o« l| 
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:i ~^ ain (X log + B c08^ (X log 71 ' (32")' 

where 

7i = 1 - ay (17) 

. -/(^)^-^ .(;^^)^ 

The "boundary conditions of- equation (14) then require 

that ■ - < 

B = 0 (34) 

and 



tan 



X log (1 - a) I = SX . ' CzsX 



Equation (35) has an infinite number of solutions which 
in conjunction with equation (33) determine the critical 
values of 3 for which a solution to* the problem exietst 
A numerical evaluation of the solution is presented for 
two degrees of taper: 



.^t ^f^y • tip stiffness^ e^l X, ^'"^f J -If'N " " ->'• 

Irf th 1 8 ^r^' -root »<itiatixm ♦ K Z&h'y ^Iritet^* i » 

heoofll»:»- > C 'iO/ 1 • -^a * •.!i^.'!v:xa5i-iT:?o jt'^lG-.;- .r: .V^^;^CiI^^i> 

is found by a conventional meth(itf ^of 'VuccVWVVe^'^ 
mations to be 

. ^} w X s 2; 546 

and equation (33)j t.hea g:ire,a .th;d>)fl^cm?,8j>.o^ding critical 
value of B, 

3 » l*a97. 

"t^^^iiiiri^ri^^ Q^l'Xifli^l^V^^^^ ^^e., wing ,^4 , v'j 



(36) 



/v^ tip atlffne»3 1 / 5N 
- ' root atl^ffwetf* ' 3'6' . . ••V^' -"iy^; v" ' 

lor this oa8 0:,t whloh ,wa9,i«T.alu^1>:«ftr 

solution for the first mode is obtained in the forn 



3 *'1.016 



(37) 



r(y) = .« -ein- ^r;ii2 log' (^-f y)] 



6^ 

1»y «(iuation (19) where u-'l/S' Sritf -ff -'iii 

functions of ordsr ±1/2 are expreasllile In terns of tMe 

pl,r*ular fnno-^lone, so that the general ieol.jitlon of ooua- 

It th« iJOTindary condition *(0) « 0 le impoBed, it^f ol- 

' V it ...... 

' *(y ) « 8ln I (iO-ny|) . ^ j tin Py ^ (99) 
Th« ^oujidaty coiftitton ♦*(1) » 0^- detexalnea tl»» «*Xti««l 
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IquatiOBs (13) aiwi (3*). dBt^tBla*- th« auxlll*?y Uf-t 
function r in the form 

r(y) « A «ln Py 



In particular, ir a m..l/2; .o Mhat tK'« tlT> chora is 
one-half the root chord and th^ tip stiffneea 1« 
fourth the root 9tiffne««, the flr.t deflection aodt of 
half the wlrg le deecrlTjed hjr lihe ea\latlon« * j 




(40) 



(41) 



P 3 2.039 
r(y ) = A «ln 2. 029 y 



i Jy 



rT.ela.2^029 y 



(42) 



1>y-a%«aUon (33) in tha for* >.v , ; c^J^i^ro 
y,^a I • 



vhera 




(44) 



while the houndary conditloae of equation (14)j«cf 
fiad if 

(45) 



asd the paraaater X le a eolutlon of tht aquation 
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In the case a = i/2. where the tip chord ie 
balf the root chord and the tip etlffnes- ^ ' 
the root etiffnesfl. the following data are ohtained for 
the firat deflection mode of half . the wing: ,,,, 



9 = l'.'653 



I 



as.'..' 



(y) = ^ eln r2.946- l»e" (x..i4 y-^^ ^-■ 



,ln [2.946 log (r- ^]^ 



d(7) 



(1 



Solution hy a Method of Succesalre Approilaatlons 

' The detetfttlnatlon -of the eBalleet critical t«I>x».:o£ 
the parameter P- and-of the oorree]oendln<f functlone^^.J,, 
7nl T froB equations (12) and in ^^^^^Z^" 

solution of equation (15) In deed form is 
oMained.l. Conveniently, accoBpllshed hy a jj. . 

.Juece.elve approrlmatlons elmllar to a method aeeeclatad 
with the names of Stodola and. Vlanello. 

A oonTenient function' O^ly) ^ % ""\;5?!;oni%' 
such a way that it Batisfles the boundary o^nditlone o- 
JquitJon^di); A function rii(y) is th.»^4t*er.ip.d by 
l»4r«du«J|i«.,thl.f t^ipctlon, in^Oj equation, ii?^^,,, ., 



I M '^ft z"^* -f.'^il a':- 



If the function 



it introduced into equation (13), 



a: i « 



(49) 



the reeultant equat_l on.V.,.;i,;;' 

flan he e»lT»d for by dl**eirinte«ration.,»c y 
ary condUlone of •qu»tton (14) 1» »y 
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If 61 (y) were the exact solution of the problem corre- ■ 
spondine to a critical value of P. it would be POJji^l* 
to chooae & so that the functions *i and «i are 
Identical. When thi. is not th. case. 

to the dealred critical value of 0 can be i'^^'^i^'/^^' 
it Is required that the function?. and Ti agree as 

well as posaible over the Interval |y| « i. This J'*"'^- 
mination ie usually accomplished by requiring that the 
two functions coincide at a suitably chosen point y = yo . 
so that a first approximation to P is given by 

^a^MTol (51) 
^ cpi(yo) 

A more accurate procedure proposed here consists 1^ 
QUirins that the integral, over the span, of the dlfj"" 
•nce between' the functions *i and Ti vanish, so that 
a first approximation to P is given ^7 ' 



S 

Pii = - 
S 



*i(y) dy 

(63) 



If a convenient multiple of *i(y). is treated 
second approximation. ^3(7). process can nov be re- 

peated indefinitely, and it can be shown that the llBli- 
ihg value 

. • '• ■■■.1- •. . • 

3 X-i «n(y) *y ■■silver- . 

P o lim ■ - 

elves tht smaXleBt critical value of vhlla th# li«it- 

Inr "^^ia^^i 0^ ^'^ ^f^'.^** - .r^*^ 



j«i^,»»«fl^f *^^ii-.'*r 

/reprtsents the corre»l><>ii41n< twi«tl»« aod» of th4 "Inff* 
l81no« th« amplltud* of Hh« function « It lnd«t«r«lnat« 
vithia the fraaeworlc of the. linear t^^'or? » f LI" 
lent at the beginning of -*aehv eyol» to Bi^hlfr^tJ^^ «**^ 
proxlmation glT»n 'by the preceding erol* ^^^^ 
the Initial approiUmtlien- ♦» in- each) dycle hae a 
Bum anplltuda of unity. ^ « 

Hapld conrergonc* of the procena hae heift-^t****^ l*"'" 
the Initial approTimatlon «j(y) le.deflned ae a eulta- 
l)le multiple of the eoj-utlon of equation (13) oorresuend- 
^^ing to a unlfora dletrlhutlon of llf%>long the «p*a, eo 
that 

/r.. Ae an lHuitrailon of thle procedure;- tHe apeclal 
'oaee of a unif orm lining for ' whlch v 

c»(y) - e'(y) - i'(y) - 1 (5i) 

t«>an»ly«ed. Ae before. It le- eufflclent to oonelder only 
one-half the wing." ihe Initial approxlaaUon to « It 
deternlned hy replacing fij) "by a conetant In equation 

illi . constant (86) 
dy» . 

n -. ■ ^ ' • 

The eolutlon of thle eqijftlon eatiefyio* eondltlone 



« 
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The oonetant ie ar^bitrary and ia conV^nlently cho«en ao 
that 1^1(1) =» !• It then follows that 



Troa .-equation V (48 ) therf» follows » 

. IiCy) - *i(7) = 2y 7' (59) 

and the introduction of Ta(y) into eqaation (49) yirlda 
the equation 

' ' • . ilii. « P^y'- ay) . (60) 

dy^ 

Integrating equat i on .( 60 ) twice and imp08in£ the 'boundary 
conditions of equation (56), the function *x deter- 
mined as 

= ^\fr • (61) 

The condition • • . ^ . . ' - 

1 1 

. ■ *^ > 

then gives the firtt approximation to the critical ta1\x4 
of 3. 

3 J . = S, 5 

> (63) 

Which differs from the eract yal^ua^. : " f ' I'^'OS, "by 
0.66 peroent. 

If th« cycle ia repeated, starting, wittt the ?.ltfitrl*I 
approxlmatloa : ' 

♦a(y) » - i- (8y - 4y' + y* ) (64) 
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it is found that 



and the condition , 'f* y 

gives a' Becond approximation to P," • i^^T^^'r*'^' 



Pi^' - 17 1 
« 1,6718 



(67) 



which differs fron the exact ralue hy less than 0,07 per-* 
cent. 

Since the initial assuaption for * in each cycle^ 
was so defined that its maximum value (at y « 1) is mxity, 
an estimate^ of the rate of convergence is afforded by 
comparing with unity the maximum values of the apprexi- 
aations obtained after successive cycles^ Thus, in 

the present example, equations (61 ) and (63) give 



xi) « £ x -^ • 1.042 ^-i*" 
and equations (66) and (67) give 

If the succsssive approximations i« P were 4'#1ri?t^ 
; mined by the conditions 



Ml) « *n<l) « ^ 
-4a place of the integral ,C;0:nAiii^ used here, the^ T^tffe 
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p. zr 1.5492 ! 
p3 = -1.5681 j 



(59) 



which differ fro^ the exact value hy 1 38 percent and 
0.17 percent, respectively . would te oljtalned. 

A further illustration of the advantage of the in- 
tP^ral condition (eauation (52)) over the condition 
u ually imposed (e,uatlon (51)) ^part 
alderation of the fourth case treated in the earlier part 
of thfs section (equation (47)). In thi. case the ini- 
tial approximation to t is of the form 



^(y) = - i (yf' - syf' * 2) 



(70) 



2 

where 

• 1 iv (71) 

and, using equation <48;) , eauat 1 on (49) heco«e« 

^ ry;* ^1 = 1 p' '(yr - 3 + 27a»r (73) 

d^i L dyi J ' 2 - 

If eouation (72) is aolved. subject * % J^^^^S" 
ditions of eouation (56). and. the result * °f-f .^^/i;3t 
equation (70), is introduced into equation ^^2) the first 
approximation to P is ohtained as 

=V.664 • " ' (73) 

This value differs from' the exact value 0 » Jl^^iJt'J? 
in eouation (47) hy ahout 0.66 percent. If. In place oi 
eouation (62), the eqUBtton. _ 

?,(1) = 1 ■ 

is used for the determination of 9. an approximation . 

l9 ohtained which differs from the exact Talue hy ahout 
2.48 percent. 
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PROOBDtJRS roa SOLTOIOH OB TE* PEOBLKM ACCOSDISO . 
TO illTISO^tlNl THlOBt 

. . . . / 

While direct treatnent of equation (11 J by approx- 
ijaktV method^ le' V^ffl'Vle, it* ir alore conr Anient ♦ part^io- 
^^IkrX/ iji dealing Vitix thr ortti.eal defleetloh mode#,"tO 

Vd¥lt vlfth nuatTont t4) jthd (-6) iw'd tc? proceeil a »e^ti- 
od or eucceVeire approxiBatlone^ ii«ira¥ «d -th^^^ :i^*otfa4tfre 
given in the preceding eectionr Ac^co^rding ttt'thef It flying- 
line theory the' tymmetrical and anti8ymaetrical Reflection 
motfir WrVeepcmd. in general,' fea diffe^reVt '^rWoefl fllgiit 

Ve^l'ocltier'." In* what folloVa , the treatmeiiV ie rireti^dat^ed 
to the symmetrical case and treatment of the antieyrfte^t^i-" 
cal' caae ie left for future work^ Attention then may l)e 

vjf#«trlcted at be.tore to, one-haXf the e^an^ (0 < y < !)• 
ll%o, at in the preceding tectiong it it aetumed that the 
wings considered are Initially «t. a f ero angl^ of attack^ 

. 8.t ari i ng 'with an ''inl t i al * is tuW* i * i (y ) ♦ * dct«x^ 
mi'iSed 'aa Vefore as the deflection cor reeponding -to d uAi- 
i^orm dfatrtbutlon of Xlft -along the span^- according tQ^; 
the differential e(iuatlon , ^ 

c:) j^i»(y> l^J.. (const,) X e*(y) (76) 

^an^.th^ >oundax:^^ v^i^O) * ^I'.U,^ " ®r 

functlCtt. 'JTjtjr) it next daftaed hy equatlpji Xs);*: ivct*"^^: 

E .Z" liL *,(y) ' (77) 

cMy) -^^^^^^^ 

and the boundary oondltioi^f i^:,^:,, 



(?8) 



i|,»l|t« aolutlon of thU ♦^uftlon ie oomrimltatly 
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and the parameters B, . . • • • ^9 are determined by a 
ceUio'd of least eouares. The flret term of the .approxl- 
natLon of equation (79) ia included since, a 3. Is shown in 
reference 5, the contri\:ution of that term to. the induced 
angle- o.f attack has- the, required di s-contlaulty. .la Its.,, 
fi-r-st- dfirivative at the root (y = 0). • •. -•. .^if, 

•. , . . • . . ■• •' / i 

•v" If the. fu-nfition T^iy) determined from equation . .. 
XT7)- -is ':iii'*roduced into aquation (4). the . resultant, equa- 
tion - - • ^ * 



_d_ 
dy 



l'(y) ^1 , - p%'(7)- ?.,(y) • <&0) 



s.ubiect to the boundary conditions of equation (8^ can be 

solTed for by . direct .int egratlon , after vhlch the ^ 

first approximation to the critical vaLne of 3 . !• 

determined an before by the condition • 

' P -^ti^y) dy = r ^li?) dy ^ (81) 



A check on the accuracy ' 6f 'fhis approxl'mation can l>e 
had by comparing ^the maxltouB valuesiof the funetloa*; «i 
and ^1. Thus if *i Is chosen in such a way that 

• •-i),-<jii Hf; (83) 

the degree of approximation attained is indicated toy ^the 
closeness of the approximatioji ' - ■< '■ • - 



(83) 



paper satisfactory results are afforded by a tfagl* eyoi 
of operations. That is, the initial approxlaation «i 



determined from equation (76) it sufficiently eimilar to 
the exact deflection mode .that if the condition of equa- 
tion (81) is satisfied^ the two functiona *i and *i 
agree closely over the entire span. 

In case only^ a ^first approximation is required, the 
value of the integral / *i(y) 4y which is needed in 

the determination of ^^in he found without explicitly 

determining the function ^i(y)» . If equation (A) i%.^.;v^^^ 
integrated twice and the "boundary conditions 

*(0) « I*(l) d'Cl) = 0 (84) 

■ • * • " : f « . . • - • 

are satisfied, the resultant expression for ^ can he 
written in the form 

*(y) = J' 0(y.^) e'(Ti) ff(ti) dti (85) 

0 

where, the Green' a^functipn.. ft . is. defined Jby. the. eqtuat ions 



(86) 



lind ^^ 
It th«n foUowe from equation (86) that 







•1 

7 


0(y,^) «= < 








i ,«(7) / "^^^ 


1 
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.1 



G(y,^) dy = / g(y) dy + g(Ti) /• dy 
=« / (1 - U g'U) d^ 



o 



^ " ^ d| (89) 



equation (88) can "be written In the form 



?"i(y) dy 3 / H(y) e'(y) ra(y) dy (90) 



where y 

H(y) = / V" dti . (91) 

If equation (90) le Introduced Into equation (81), 
the relationship determining the first approximation to 
3 can be written in the form 

P ' „ 2 (92) 

1 

H(y) e*(y) ^^(y) dy 

S^uatlon (83), which afford* a check on the accuracy 
obtained, can also he expressed in terms of the function 
I]^ if it is noticed that, from equation (85), 

*i(l) « Pa'/ 0(l»Ti) o*(n} F^(n) dr) 



symmetrical deflection mode can^'bu'^iimmari'se'd' as" ffSi'iwin 

(1) Determine from aquation (76) and the 

boundary eondlt lone* t ^i-(^)' » ! ^-iX) V>::0^v ^nd. deter- 

mine the multiplicative const unVjto $hat ^x(l) « 0. ^ 

, (3)r Petermin#7 (y) ^. ^>ir-0|i ,eq-gtat i on (7?) and the 
houndary' oonditioni» . TiCil) « Q hy the procedure of 

(3) Determine >f rdir' '»qua^ i^n (92h A d^^^4c^ otf - 

the aocuracy of the determination in provided by equation 

(cWy< If 'greater accuracy is deelred. the function *i(y) 
le next determined from eouatio]^ (80)^ and :.«the bj^undary 
conditions h^iiof « t^l) ' (]a^'= 0; A ,fiinctlbn da(y.) 

is then defined as 

^^V/ d^(l) r itnrio-;- ^Mti- 

• 0 that ^|jfcM^>-?'A>rU>$^ ;«t?^cf;^*J^i?*>^-: '^^'^ 
determined as in step (3) of the preceding paragraph, and 
a second approximation to P is determined as in step 

proad^^fc;-^^^! ^ft t.h.^.44i6lri^lT*fivP^S«i^«^ 

tion d.x^-r:.JBM^r^ c :.;.:v^..c,. :vO/ ^SfcX^ltf Jhr<HSbMM?»^3^ : 

at nine ei^iaiyepacei point t.alQng the sealepaa« The 



26 



NACA Technical Note'lTo. 926 



▼ aluei of the function at the same nine points are^ 

then determined from the data presented in reference 5 by 
a purely numerical procedure requiring less than two hours 
of time. If the values of. the functions e*. H, and g 
are tabulated at these points, the integrals needed in 
equations (92) and (94) can he readily evaluated by 
Simpson's rule , modified if, necessary so as to take into 
•account the fact that the function l'i(y) has an inifi- 
nite derivative at the wing tip. (y,*. l). the approx- 

imation of equation (.79). is used, such: a modified formula, 
derived in the appendix/, is of the form. 



v(y) yi(y) ^7 



) ^'k V 

"+ .0..00506 -Tb' + ^.Ajjjj: (95) 
n«o 

where, yjc ' k/8 and a.^ is the Simp e on ' 8-rul«- Weighting 
coefficient asaoclated with the. point yjj 



? 24 

• r ■ 
9B 



24 



24 



24 



24 



2_ 
24 



(96) 



The last term of equation t95) gives the correction due 
to the fact that ya(y) has aA infinite derivative at 
the point y » le 



LITTIHd-LINS A^AltSIS CFJ KCPLICIT' CTASia" 

In order to illustrate the procedure of the prec^'«^ 
Ing section, and to investigate the importance of the 
aarodynamio ipati- effect , th^j ^as^^e «bn»ly^ried tnv«^ 
Ing tedtloft according to the section-forcflif* thaO!»» -W^ri^ 
reccHildered oh the baeie of th^i llf ting^^ln^-Wiilwy •> 
«lno#' oaXy aymmtitrical deflections *r© to^ tTem%«a^'5 
etrt^^ton^iillt be r#ftria.t^d t».^ii«>-hal4r :^e^'«rlngi^,v .noic^ 




In the numerical calculations it is assumed that 

4 h 4 

For an untapered wing this value of n corresponds to an 
aspect ratio .-..^ Ti :»»nv: i^ u ;^ •: 

1 = m ( « 6) - ::;c-- ■(98) 

c 

while for the general case the corresponding aspect ratio 

h /S. where S is the projected wing area, is g^v^a "^7 . 



h 



a 



(99) 



S c*(y) d7 



.In particular, for a symmetrical wing with linearly 
tapering chord, 



it follows that 



L. Br— (101)^. 

3 1 - i a 



It should "be remarked' that the calculations contained 
in this section were made hy retaining a larger number of 
significant figures than are incliided in the data presented, 
Thus if the caliulationa are repeated ok the. ^asis of the 
tabulated data, the results may differ slightly -from the 

The auxiliary functions H and g are determined from ^ 
equation* (91) and (87) in the form 
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H(y) = y - iy» (103) 

and 

g(y) = y (103) 



Ab in a preceding section, the initial apDrorimat i on 
to the first symmetrical mode (for 0^ y ^ 1) is deter- 
mined from eauation (76) as 

^i(y) = 2y - y^ " (104) 

Next the function is to "be determined from equation 

(77) which hecomes 

1 

J^(^y)^l.f i!l , 2y . (105) 

^'^U^l dri y-Ti 

If the approximation of equation (79) is assumed, the pro- 
cedure of reference 5 determines the constants as follows: 



B a 0.4298 Aq » 0,2680 



A^ » -2.3442 A^ 

4 O 




(106) 



The data needed for the computation of 3^ according to 
equation (93) and for the check calculation of equation 
(94) are listed in table 1(a). Thus, according to equa- 
tion (92), the quantity is the ratio of two integrals 
of vhich the first is obtained by Simpson's rule as the 
weighted sum of entries in the fourth column: 

/ *i(y) Ay ^ y *i<yic) « 0.6§67 (107) 

and the s^con^ oo^:a^ned in rirtue of equation (96) ae 
t.li-> ::.^?» c.f %*wij:tited ^xwc^^cts of corresponding entries in 
tlxQ set?, ad .'^ w riVih coluane, plus a correction tarns 
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(0,00506)(1.4773)(0.5000) » 0,1654 ' (108)- 
Uquatlon (92) then glv«ii : ' • • '^^^^^^ ^:t.:at^-£> 



» r u.^..:.i r y.:^ 1 . .(109.) 



In the same way, the i'ntiegral ^contained 'in equation (94) 
is evaluated as the sum of weighted products of corre- 
sponding entries in the. third and fifth columns of tahle 
I (a) , plus a correption teriji 

1 ^' 

4- (0.00506)(1,4773)(1,0000) 0.2479 (110) 



Bquatlon (94) then gives 



(Vhlle In the present cai*--ilie^' integrals in equation* 
itOt), (108), and (llO) can ea«lly "be evaluated directly, 
the numerical nethod of. bralttAtl«n ^^uat outlined I* par- 
ticularly convenient In Case the Integrands hare compll- 
catod analytical expreB»is?;i_^,,pr: ,aTe.,^^^^^ ' 

Since equation (93) reoulret that the Integral of 
ih* -difference between-tl8t*£iwea*«i»ST* »PP'f**^"^ ** 
and . «^ \>a Eero and equation (111) «how«> that the two 

elnce hp^h.,fi»»ft1'l '«'«W5^f^^ T».^t.vaad,h.aT^^|^5ft««jj<i>^ 
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• lope at the tip. it would appear that the agreement T)e- 
tween tho two functions over the span ie euch.that the 
process need not be continued. The" functions and 
r, would then he considered as the deflection and lift 
modes corresponding to the critical value given in 

equation (109), 

In order to verify the accuracy of this approxima- 
tion the second approximation is now determined.' The 
function ^i(y) is obtained from equation (80) hy direct 
integration and is tabulated in the sixth column of table 
1(a). If the values of this function are divided by tho 
value of the function at the tip. the corresponding values 
of the function *a(y) » 4i(y)/<i(l) are obtained and 
listed in table l(b). The parameters specifying the 
function ra(y) which satisfies the equation 

^il. *3(y) (112) 

dn y-'H 

'\ 

are then found by the procedure of reference 6 as fol- 
lows: 

B « 0.3884 Aq » 0,2624 = 1.4074 

(113) 

A^ « -2,4319 A^ « 1.8538 

The values of the function are pretented in table 

1(b). Jrom equation (92) the eeoond approximation to P 
is found 

and equation (94) gives the reeult 

ii^ tir. 

Tor f ixi»l ooaparieoir, tit# :?T«taw ?fiaWttW|;C^4[Cj4«o 
deterirltted: from ' Ta(y) £t i^^itWt^ii-im 
equation (80)^ are Included Itt %%^X% !(>)• 
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' may be f .lJ^*J^i;C?^?^T*T^ 
rapid rato of < ooaYarsettoe, .a|,. the ^J^t#2f,at4jy a. jirocet t ocours 
if, instead of taking aV tire-initial awMPOximation-to thi- 
de f 1 e 0 1 i on m oitf -:t h e e pLu^i , o f :^tqu«tti.d4 76); « o lu-- 7 
tion of tha proxies apoording to the edotionrforoe theory 
ia taken as the initial approximation;* The r^ft suits of a 
oaloulation hased on this propedt^re,. with*^^^ ^ , , 

x'i *i(y) « *,^(yJ « sin | y. ^ vU?^ 

are presented in tahle 1(c)/ While' the >alue obtained ' 

tor 



7 



« 2.004 (117) 

and the values of the function *a(y) agree cloeely with 

the preceding results, appreciable differences are pres- 
,ent between the successive approximations *i and Zi 

in this procedure, , An^ indication of the fact that the. 
function giren in equation (116) does not afford an ac-; ' 
curate 'approximation to the actual critical deflection^ 
mode would be afforded, without a complete explicit eral*- 

uation of the function by the readily calculated 

Talue.of djd) for this solution, 

=r 0.9682 (118) 

yroa thase results three useful conclusions , which 
will be further substantiated in the follewing treatment^' 
may be drawn: 

(i£X5 The approximate deflection mode determined as 

the solution of e<^uatiQn (76) is more nearly in agreement ;3 :^ 



with tie actual mode than is the m ode predicted by the 



* y ?! ' ; ./ 



sect i 0 a-f or c e .ilbe ory . * ' ' " • • * ' ^ "'^ 



(# %he preVe^t pithbr hlfj^l^ -'^^ 
s ensitive to the choioa of tha initial approximation io . - ^ 

thVli^i tl ^Oti^; ko^e as''Tt>g^^ ^ the 

f ir i^feoxVisati on^ ^6 thV oti%mtVinm of ^"y; ""^X'-^ 
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(3) if^ gj l« determined by aoustlon (9 2) th<- degre» 
of approximation attained In equation O*)' &ffOTi.n a 
reasonably accurate estimate of the agreement betwea the 
Initial approxlteatlon and the actual deflec tion «ode. 

If the reeulte preaented in equations (31) and (114) 
are compared.. it is seen that the aerodynamic span effect 
il responsibie for an increase of about 28 P«"«f 
predicted ralue of the torsional-divereence Telocity. In 
figure 1 there is presented a comparison f 
Ipp^oiimations and ^3 to the deflection mode, 

and of the sucoessiTe approximations and I, to the 

corresponding lif t-diatrlbutlon function. The mode 
*,f = ^8f - i y predicted by the section-force theory 
is also included. 

2. Uniform chord, quadrat ically dec reasing stiffnees 

[c ' ■ e* = 1. " (1 - ay)']- -" The functions «i. H. 
and c are obtained from equations (76), (9l).. and (87) 
in the form 

i- . u - (i^ - 1) m.) 



H(y) - -^r^^e ,^ - 




(130) 



«her« 71 - 1 - ay. These function*, to^rether with 

function- rx(y) determined •^^•^i'*"'^"!*.?^* 
uat.d for ; - 1/3. in table aCa),.!**,,^©^,^^ ^d f^^ - 
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^^J-yi J^7M9i:'^-*3 &i''>^' •v?'?*;^.'-®.]"- ^*«■VW:7.5.«• 



P « JL 708 I 

^/r^ v^5> '^^'^c .}oi ; ^- • vj (SS^/ ^•^^'^ ('il^v 

Irom equation (94) there f ollova f or a^« 'l/3, 

-I't^iift A.-Ti?JiiC?T3'3r -;;;;it' ^ >^•^.^'^x^/ ^ 

.... — T^(iy:iriij9m'^'^ 77~ 7" \ (X3^j|| 

and, -ford si:«i:5/6-^ ^ - ^ a .^>t>^*^q -rii^i :r ^rK^ i tvr.r^^t 

7,(1) » 0;^8»7-'^ - :r^o;o- <f25)^ 

, > ; ^ ^ ^ 

SinSe*'equa!ji^nt ri34)'and (125) indicate iTeatlB^f act ory 
agreement "between the functions and A, in "teoth 

casest it is concluded that the functions^ . and., 'i.. ^-^v 
af ford. ieaeoa«(>»le 'approxiaatl one to . t.he .Qj:.iti|aCal 4*'-^.fA"' ^ 3 ' 
tion anddff t ' modes' and that the comput ed ^Yalues^ 'of 3 
arer^Wf#ia^t#^«lfy ^i^cruiiat'd/^^ ^^^^ .^'-^ ty^*^ 

->-JiXf^ t«n'Ee^^tf«tf«e^'^a"^^f/3. the* ixii'^lAf ^V#4xli^^^ 
ie taken to 'be proportional to the def lecti on^odV ^'''^ 

predicted ty the section-f ores theory (equation (36)), 
tl^tjttalue 

ie obtained and ie seen to he in good agreement vith the 

>glTen in equation Cl32l)^ ^Howja.v^. the check cal- 
cutatrion in tbit case gives Vhe resVli; 
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which result , if compared with equati on (124) , aigain in-'^ 
dicatea that the deflection mode doee not afford so. 

accurate an approximation to the true deflection mode as 
does the function defined "by equation (76), 

The functions *i and are represented graphi- 

cally in comparison with the function ^gf f ^gf P^®" :. 

dieted hy the sect i on-force theory, for the cases a = l/2 
and a 5/6,' in figures 2(a) and 2('b). If equations 
(122) and (123) are compared with equations (36) and (37), 
it is seen that the aerodynamic span effect is reaponsihle 
for increases of ahout 32 and 39 percent in the values of 
the t orsional-divergence velocity for the cases a » l/2 
and a 5/6, respectively. 

3. Linear chord, quadratieally decreasing stiffness 



al~ay, (l-ay) 



3} 



The functions 



and g in this case are identical with the corresponding 
functions in the preceding case and are given in equa- 
tions (120) and (l2l). From equation (76) the function 
is determined in the form 

^i(y) = ^ ^(1. - yj - (1 - a)^ " ;^)} ■ ^^^^^ 



where y^ 1 - ay. The values of the functions e*, H,^ 
g, and as well as the values of f ^ ' - dcit erilried 

from equation (77), are listed at tha nine po^n^t^a ^aa<^d«4- : 
for the approximate integration in tahle 3 for the case 
a » 1/2. - Jlquatlan J92) then determi^iei the -flret approx- 
imation to .fi, .. . ^ . _^ 



Pi' ^ 5.637^ 
01 « 2.374 J 



(129) 



while equation (94) gives 

■ •■T.d) ..0.9953, J 

Equation (130) indicates that the reeults of the first 
approximation are sufflcientl-y accura.ta to juatify tar- 
alnating the process. 



ta* •H.lt infei-lii* tie ory. pt^dl a^ft :A?,d i'Wjarg.wiOi*. ▼ •l oci ty 
Ctht* a*«* WHtcH 18 about; a^cp«rcetflt .>4«fe».r thaJk>l»«- 

tlon. aa4 lift »od«» acftoaAii^ i^o ^tii*, tfe««^i:*%v^%(«e'i* 
compared in figure 3. .■ ■ ' > ;\. > r/ -'>:.':'••• : t--i^.*::- ^r ^* J "».'. '^ 



6a 



vhoro y'l " 1 - ay. The functtone needed for the coaputa- 
tio» of 3 are oraluatcd in ta\)le 4 for the case a « l/2, 
Trom eauation (92) th;. appr6xlJiaatio;n 




(134) 
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Tilt lift and deflection aodee *i: and. art cob- 
pared with tha corresponding results given by the. eaction- 
foroa theory in figure 4. A comparison of equatione 1134-; 
and (47) showb that the aerodynamic span effect la ra- . 
gponsihlo for an increase of 20 percent in the ralua o£ 
the t oraional-divergence velocity. 

lor thepurpo«e of further verifying the accuracy of 
the present procedure the same case has heen analysed hy 
two other methods. lirst. if the initially assumed de- 
flection mode is taken as a suitable multiple of the mode 
predicted hy the section-force theory (equation (47)) the 
approximation to 3 is obtained as 



=» 3.898 



« I.974J 



(136) 



while equation (94) gives 

did) « 0.9332 



(137) 



Second, if the assumed deflection mode is taken to he the 
deflection corresponding to a uniform distribution of 
twisting mbmept along the span (el » const-^) rather than 
the deflection corresponding to a uniform distribution of 
lift (I « const.) it is found that 

a 1.976 J 

and 

« 0.9567 



(138) 
(139) 



▲ comparison of these results with equations (134) and , 
(136) indicates that the deflection cot^responding to d 
uniform lift distribution is in closer agreement with fhr 
actual divergence mode than is the deflection correspond- 
ing to the other assumptions. The remarkable agreement 
between the computed values of 0 shows again that Inso- 
far ae t)xe determination of (j le, cox^cerne^ ^the^^ ?r#>ent 
p^cd#dure ie not extremely eirn#ltlVe. the^^fitt4 a tW 
InitlUly aseuaed. dtfXtotlon mode. 



5. Linear chord (c - - 1 -^'.y^J 4tiW«H^>^o^ 
. .,.>i>ion ,tlffn„« (fig^.,&Kc>«.:^;A**^ application of 
the method, of thi. paper a e7»««trlcal wins i» 

in whic^ the, cl^ord .**W,- •^'^•^^^4^^^^^ Hv- " 

haLf th* roor TV««: "a:th*' sectrott /fn^^ 

:t:ironrr.i5 ::^d*ra/)Tar rc^o^i-s^d jrsr4%- 

section t,aper and in each ca,e counting ^J^; 
tained h'etween the corraeponding J'"" 5. 

If e^uatipn (76), which h^re take^ the f or« ^ ^ ^ 

is integrated twloe and the- eondlMon. ^(0) « V^^^ ' ° 
are i«po.ed. the function can De written in the for. 



..\ ;vJ. , ■■ ■■• : ■■ ■ ^ ■ y 



Thus the function can ^. ccareniently evaluated fro. 

- y + 1 y»i' 

the function ^' ^^^'^^K ^ 5'* 

t.i-iiVlUft: Jv isi.r*i»d^ .v; ^^^^ 

-.^,,«*c4ofi (77) y>T th# procedure or reference o» •t'^*^ Iw^ 

a*tl>'o*^^ir ai^I>roxHkat«. tn.t«*r*^il^tt ^'VJf*: ^f:?^;^ ' 

:-WU6?M^5.i4,#- ■^t.^.-c; - ,ri•; 



> 
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ZquatioYi (94) then givee 

5"i(l) 3 0.9979 (143) 

so that, while extreme accuracy should prohahly not he 
expected as regards the lift and deflection modes in the 
neighhorhood of the discontinuities, a satisfactory agree- 
ment between the functions and is indicated. 

According to the section-force theory the first ap- 
proximation to the lift function F(y) is given hy 



Tg^(y) s { 1 - JL y) *,(y) (144) 

in which case equation (92) glvea, as a flrat approxinsa- 
t i on , 



Pgj= » 1.364 



'1 



^ (145) 



« 1.168J 

Comparison of equations (142) and (145) shows that 
the aerodynamic span effcci) is responsihle in this case 
for an increase of about 19 percent in the predicted val- 
ue of the divergence velocity, The lift and deflection 
modes for the two theories are compared in figure 5» 



COirCLUSlON 



The results of this piaiper indicate that neglect of 
ithe aerodyaaflio span , effect may lead t,o^ cn appjeciahle • 
underestimation of the torsional -divergeccs velocity the 
d^lffeTenqe: hetween th^ values, obtained^ with, and without 
ne^iect of thj-e effect amounting to^ 1? to 40 percent in 
thV numerical examples .presented. ^^^- '^ /i^ 

^'^1 ^ IjtCVteW -^of the iiQt t^at' ttMs;Vae^ 

oialy Vings with a span or a1>out' six t^lmWe the r dot- chord, 
it seems desirable to consider' a grekte'r rikrl^ty 'of vlhgs 
and, in particular, to l^nvesti^ate the. relationship he- 
^ tw^en the relative ma^nitvtfe ^6f the aerodynamic spaa ef- 
fect and the aagnitude^^of^ t)it aepaat ratio. 
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An extension of the p-resent pro-cedu.re to- the analysis 
of antlcymmetrlcal deflection modes, aa well as to the 
analysis of the elastic deformation of wings before tor- 
eion'ffl divergeixce occurs, can her accomplished withatit- es- 
sential difftoiaty,. 

Massachusetts Institute of Technology^: . ^ 'i , 

Camhridge, iiasa.. Feb. 1943. 



APPBITDIX 

A JORHULA TOR APPHOXIMATS IS'arBGEATIOK 
If the approximate Talue of the integral 

:) dx (1) 



is required, and if the function f(x) is of the form 

f(x) = p(x) Jl - x^ (2) 

where pCitl finite at x » 1, conventional formulae" 
such as Simpton's rule fail to give accurate results due 
to the fact that f(x) has an infinite slope at x = 
A modification of Simpson's rule which takes this fact 
iito account- ftfliere derived for a nine-point weighting, 
systea* ' ' ^ '. 

. , .1fi,tjtt^tka, .notatt<)|a,^ f(k/8), Simp«oa»8 rule gire^ 

/ f (x ) dx 0 4f , ♦ 2f a 4f 3 4^ 2f^ ^ *f ^^^t^ j ' 43> 

If^ In the raa«« 5/4 < « < It the function fix) ii 
approxi«at«4 by th» sxpressl^n ! 
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fix) - ai -/ 1 - X f a, (1 - x) ^ a3 (1 - xV 



(4) 



the oonstanto aj , a^ , and can be determined 80 that 

eauation (4) i« a true equality at the pointa x = ?/4 
and X 3 7/8 and so that the derivative of the difference 
between the two aides of that equation is finite at x = li 
It then follows that 

aa =/2 p(l) 
aa - 4fs 4. IGf^ - (8 - 2 J2) p(l) 
^3 ^ 32fg - 64f7 (32 - 16 ?(l) 
With the approximation of equation (4) there follows 



1 



f(x)dx ^ a, ^ ^ aa ^ Y92 a3 



32 



or 



>4 



f(x)dx -^[fs ^ - 2) p(l)j : 



Equations (3) and (7) can then he combined to give 



f(x)dx ^ Y^<^^ ^^^k^*:^ 



k«o 



f (x) 



'x»| 



(5) 



(6) 



(7) 



(a) 



where xy^ * lc/8 and Is the weighting coefficient 

associated with the point x^ by Simpson* e rale. The 

last term in equation (8) is a correction term which takes 
into account the fact that f(x) has an Infinite deriva- 
tive at the point x • 1, 
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TAB 12 DATA lOR UITITORM WIKO 




(aj Tirat, Approximati on 



7 


H 


S 


*1 




^1 


0 


0 


0 


0 


C.2680 


0 


.125 


.1172 


.1350 


.2344 


.3040 


.2228 


. 250 


.2188 


.2500 


.4375 


.3848 


.4263 


.37E 


.3047 


.3750 


.6094 


.4779 


. 6054 


.600 


.3750 


.6000 


.7500 


.5518 


.7646 


.625 


.4297 




.8694 


.5857 


.8691 


.750 


.4663 




.9375 


.6760 


.9471 


.875 


.4Sc;2 




.9844 


.5157 


.9889 


i 1.000 


.5000 


1 1 .COO'- 


1.0000 


0 


.9994 






: 1.4773 





'.b) S0oond Approxiaation 



— I 

y ! 


*3 






0 ! 


0 


0.3624 


0 


.125 


.2229 


.2984 


1 .3224 


.350 


.4266 


,3806 


1 .4368 


.375 


.6067 


.4764 


1 .6063 


.500 


.7649 


.5533 


.7549 


.635 


.8696 


.5896 


.8699 


.750 


1 .9476 


.5803 


.9481 


.875 


.9895 


.5186 


.9901 


1.000 


1.0000 


0 


1.0006 




B + lAan 


« 1.4802 





(c) Tirtt Approximation with *a • *8f 



7 


•> I 






*3 


0 


0 


0.2399 


0 


0 


.126 


.1961 


.3730 


.3109 


.3201 


.260 


.3837 


.3507 


.4044 


.4220 


.376 


.6666 


.4435 


.6768 


.6009 


.600 


.7071 


.8211 


.7194 


.7608- 


.625 


.8315 


.5622 


.8306 


.866« 


.760 


.9239 


.6614 


.9067 ^ 


< .9463 


.876 • 


.9808 


.6083 


.9478 


.9893 


1.000 


1.0000 


0 


.9683 


1.0000 






• 1.4606 





I 






■/_H 




J^s2. 








0, 


0 


0' 


0 


0.2090 






: .125 


r.l248 


;i333 


' .1616' 


.2371S 


: ,1217 \ 




: .250 


'■. .2484 


.28'57 


.3216 


.3069 


.2569 




.375 


.3690 


.4616 


.4776 


.3925 


.4082 




.500 


.4841 


.6667 


.6265 


.4685 


.5584 




.625 


.5897 


.9091 


' .7633 


.5163 


.7089 




, .750 


.6800 


i;20oa 


■ .8801 


.6290 


.8484 




.876 


.7469 


1.5556 


.9665' 


.4910 


.9554 




1.000 


.7736 


2,0000 


1.0000 




1.0000 










= 1.4287 





(b) I = Ih(i 





' ff 


■ 


— 'ti — 






* • i 

.' ■ ■ 


f-O 

' V". • • 


i2788 
^. V43Q5 
:;'604^ 

U:^'79e6 
:i;$;-0124 
';'r'.JU;'3346 

;:r:t;-^8Cii 


■ ;4-3158:..' 

;h-i«4B6~;:; 
,;:<;b671-- 
u'-.i.'^o-43;; 

-3VSS08 
' 6.opqo ; ! 


^.0946 

iia'ao 
" ;3ii3- 

' ♦;43«3; 
• ,5786 . 
;7836.; 


0.1428 

■ ,1615^ 

iaixff 
;376r 

- ^3428 
i39r4' 

^^z^t -. 
U342r- 
•fcl 


; x^lHa?^;: - 

•; r,-:;:3i»:, 
^ ;r'.:^33a7;.. 
! :''.i47ao " 

■ . ,.64*4 

:-\"8^4i-l 
i ilvOOOO 


/ 
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• ,TAB1S,.3.- DATA TOR LINZARLY TAPSRI50 WIUO WITH 



[^0 = Cr ( 


-1' 






a 

- 1 
S 3 


-] 


7 


e* 


E 


g 










0 

.125 
.250 

.375 
.500 
.625 
.750 

.275 
1.000 


1.0000 

.9375 
.S75O 
.8125 
.7500 

.6875 
.6250 

.5625 
.5000 


0 

.12lM5- 
.2l|«l+- 
.3690 

.5S97 

.6800 

.7^9 
.7726 


0 

.1333 

.2857 

.U615 
.6667 
.9091 

1.2000 
1.5556 
2.0000 . 


0 

.1833 
.3571 
.5192 
.6667 

.7955 
.9000 

.9722 
1.0000 


0.1868 
.21U5 
.2783 
O'*+90 
.3980 
i4l09 
.3973 
.3715 

0- 


0 

.1492 

009^ 

.^731 
.6312 
i 771*0 
.8908 
;970H 
1.0000 


0 

.1399 
.2708 
.38UH- 

.5+731+ 
.5321 
.5568 

.55+58 
.5000 










= 1.2058 





TA3L3 4.- DATA FOR LISSAHLY TAPSRING VlUa WITH 



qUARTICALLY TAP3RIHG STIJTNSSS 



[ 


» - «R ( 




Ir 0 ■ 




S 3 


■J 


7 . 




H 


s 


T 




,- 




0 

.125: 
.250 

.375/ 
.500, 

.625 : 

.750 : 

.«75 ^ 
1.000 ' 


1.0000 

.9375 
' .8750 

- .7500 

,:.6«75 

' .6250 
.^625 
.5000 


0 

.1333 

.2838; 

1^0 
.8U65' 
1.0560' 
1.21+19 
1.J333. 


0 

.iteU.' 

.3285". 

.57612 

.9136: 
1.3a*»9 
2.06l<d 
3.0791 
5*.666r' 


0 

.1173 
.2U1+3 
.3810 

.675^ 
.8208 

.9l«9 
1.0000 


O.1I409 
.1567 
.2056 

. ^2691 
.3217 
.3*W!5 
i3550 
.3%i 

.o_i ■„.;._.. 


0 ; 
.0826 
. ..1858- 

; .1066 
.••.U57? 

- -.7891! 
. .9332 


O.r r ' 

.^1?26 
.21*91 

:*3^3^ 

m 

^52«9 

...♦5000. 




• 


4 " 






-1.1686 
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TABLE 5.- DATA TOR LINEAELY TAPBEIHO WINO WITH 



DISCONTINUOUS STIJTNBSS VARIATION 







P (l 




- i 

' S 3 






y 


e* 


H 


€ 




Fa 




0 

.135 
.250 
.375 
.500 
.625 
.750 
.875 
1.000 


1.0000 
.9375 
.8750 
.8125 
.7600 
. 6875 
.6250 
.5625 
.5000 


0 

.1172 
.2188 
1.0781 
1.7813 
1.8-435 
1.9041 
1.9616 
2.0038 


0 

.1250 
.2500 
1.5000 
2.7500 
2.8927 
3.0891 
3.4061 
4.3267 


0 

.0656 
.1200 
.5606 
.9043 
.9333 
.9599 
.9836 
1.0000 


0.1455 
.1582 
. 2364 
.3523 
.4463 
.4731 
.4352 
.3770 

0 


0 

.0615 
.1050 
.4555 
.6782 
.6417 
.5999 
.5533 
.6000 








B ^ ZA311 


= 1.2400 
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